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ABSTRACT 

We present subarcsecond-resolution mid-infrared images of the debris disk surrounding the 230 Myr- 
old A star £ Lep. Our data obtained with T-ReCS at Gemini South show the source to be unresolved 
at 10.4 /im but clearly extended at 18.3 /im. Quadratic subtraction of the PSF profile from that 
of £ Lep implies a characteristic radius for the dust disk of 3 AU, which is comparable in size to 
our solar system's asteroid belt. Simple models suggest that the 18 /im flux is well approximated 
by two contiguous annuli of mid-infrared-emitting dust from 2-4 and 4-8 AU with a 3:1 flux ratio 
for the annuli, respectively. We consider two scenarios for the collisions that must be resupplying 
the dust population: (1) continuous "steady state" grinding of planetesimals, and (2) an isolated 
cataclysmic collision. We determine that radiation pressure and subsequent collisions are the dominant 
determinants of the disk morphology in either case, and that Poynting-Robertson drag is comparatively 
insignificant. 
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1. INTRODUCTION 

Excess infrared flux associated with main sequence 
stars usually indicates the presence of thermally emitting 
dust in a circumstellar disk. When the stellar age is long 
compared to the timescales for removal of pristine (i.e., 
primordial) material by Poynting-Robertson (P-R) drag 
or radiation pressure, the observed dust must be pro- 
duced by collisional cascades of lar ge particles and plan- 
etesimals or by com etary deposits (|Backman et al.lfl992t 
iWvatt et al.l[l9 99). Such disks are referred to as debris 
disks. The A-type main-sequence star £ Lep (HD 38678) 
was identified as a debris disk candidate by the IRAS dis - 
covery of its infrared excess (jAumann & Probst! fl99l . 
It is one of several stars where the disk particles are 
significantly warmer than is typical for deb ris disks 
(|Cotefc Waters! fl987t lA umann & Probstlfl99l . imply- 
ing that the particles are relatively close (<10 AU) to 
the star. For example, the ^320 K mid-infrared color 
temperature of £ Lep (Chen & Jura 2001; this work) con- 
trasts with archetypes such as the disk of (3 Pic, which 
has a mid-infrared col or temperature of ~ 160 K for dust 
80 AU from the star (|Telesco et alJ ll988). Furthermore, 
£ Lep is one of the few such stars with warm disks that 
are close enough to us (21.5 pc, Perryman et al. 1997) 
that ground-based observers can reasonably expect to 
spatially resolve the disk. Two recent observations of 
the disk drew our attention to the £ Lep disk and mo- 
tivated our program of mid-infrared imaging at Gemini: 
a strong constraint on the size of the disk at 18 /zm at 
Keck (Chen & Jura 2001, hereafter CJ01) and a Spitzer 
MIPS photometric survey of A-type stars showing that 
the infrared excess of £ Lep is high among stars of com- 
parable age (|Rieke et al.l 120051 ) . CJ01 have noted that 
the disk radius must be smaller than 9 AU, and there- 

1 Observations were obtained at the Gemini Observatory, oper- 
ated by AURA, Inc., under agreement with the NSF on behalf of 
the Gemini partnership: NSF (US), PPARC (UK), NRC (Canada), 
CONICYT (Chile), ARC (Australia), CNPq (Brazil), and CON- 
ICET (Argentina). 



fore probably comparable in size to our own asteroid belt, 
which makes it unique among debris disks discovered so 
far. As we report here, we have imaged £ Lep with Gem- 
ini to constrain more tightly the disk size and further 
examine dust production mechanisms. 

2. OBSERVATIONS 

We obtained mid-infrared images of £ Lep (HD 38678) 
on 2005 February 3 (UT) at Gemini South (Pro- 
gram ID GS-2005A-Q-2) with T-ReCS (Thermal Re- 
gion Camera and Spectrograph), using the broadband 
N (A c = 10.36 /im, AA = 5.27 fiia) and narrowband 
Qa (A c = 18.30 fim, AA = 1.51 /zm) filters. The to- 
tal field of view is 29" x 22". The data were reduced 
with the Gemini IRAF package. The flux standards 
e Lep and a CMi were observed before and after, re- 
spectively, the group of target and PSF observations. 
Flux measurements of £ Lep at N (2.03 ± 0.19 Jy) and 
Qa (0.96 ± 0.06 Jy) are consistent with previou s results 
(jFaiardo-Acosta et al.lH998l : lChen fc Jurall200lD . We ob- 
served a point-spread-function (PSF) comparison star, 
HD 42042, both before and after the target observations 
in each filter. HD 42042, an M star located 6.6° away 
from £ Lep, is ~10 times brighter than £ Lep at 10.4 /im 
and ~6 times brighter at 18.3 /im. We divided the PSF 
and target integrations into pairs of nodsets to examine 
the potential effect of airmass on the full-width at half 
maximum (FWHM) value, and there is no direct corre- 
lation. 

3. RESULTS 

3.1. Source Size & Model 

In the final stacked images, we find that the az- 
imuthally averaged profiles of £ Lep and the PSF star 
are identical at 10.4 /zm, whereas the source is clearly ex- 
tended at 18.3 /jm (Fig. 1). We have partly characterized 
this extension by comparing the full width half-maximum 
(FWHM) of the profiles determined from a Moffat fit 
generated in IRAF. At 10.4 /zm, the fits to the PSF star 
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and £ Lep profiles measure 0.31" ± 0.02". At 18.3 /tni 
the profile fits give 0.54" ± 0.02" and 0.60" ± 0.02", re- 
spectively. Even though we consider the FWHM radius 
of the profiles, we have determined those values from an- 
alytic fits to the full profiles. For £ Lep, the indicated 
uncertainty is the standard deviation of the mean of the 
series of FWHM measurements made throughout the in- 
tegration set, as described in greater detail below. For 
the PSF star, we have only six FWHM measurements, 
and the uncertainty given is half of the excursion be- 
tween the maximum and minimum values. Like CJ01, 
we do not resolve £ Lep at 10.4 /tm, but we do resolve 
the source at 18.3 /mr (Fig. 1). 

To test the robustness of the latter result, we examined 
the FWHM of the PSF star and C Lep throughout the 
integration sequence. Pupil rotation, incorrect guiding 
correction, and changes in the quality of seeing during 
long observations on a program object can result in a fi- 
nal image degraded by lower-frequency components that 
are not accurately represented in the PSF determined 
from generally shorter integration times. Thus, we have 
divided the total 18.3 /im integration sequence into bins 
of 120 s of on-source integration time and have measured 
the width of the Moffat profile fit for each of these subdi- 
vided images. The displacement between the mean pro- 
file width of the PSF and the mean profile width of the 
source implies that the source is resolved. Quadratic sub- 
traction of the derived FWHM value of the PSF star from 
that of £ Lep, as determined from the 18.3 /tm images, 
indicates that the average disk radius is 0.14" ± 0.02" , or 
3.0 ± 0.3 AU. CJ01 estimate a best-fit value of 2.15 AU, 
but they note that the uncertainties in their observations 
allow for values as large as 9 AU (CJ01; C. Chen, priv. 
comm.). Our observations made under nearly diffraction- 
limited conditions imply that the size of the emitting re- 
gion is indeed comparable to that of the asteroid belt in 
our solar system. 

To constrain the geometry of the £ Lep disk, we gen- 
erated 2-D models of the expected 18.3 /mi brightness 
distribution with a range of disk parameters. Subtrac- 
tion of the PSF image normalized to the peak brightness 
from the image of £ Lep shows a residual ring of emission 
that, given the asymmetries evident in the PSF itself, is 
consistent with the disk being approximately azimuthally 
symmetric. We therefore assume a face-on disk geometry 
for our models and calculations, but this assumption has 
no effect on our conclusions. (Our observations provide 
no constraint on the disk orientation.) 

Our models consist of a central source (the star) with 
either one or two annuli of dust emission each with uni- 
form surface brightness, where the constant total flux 
determined by observations at 18.3 /ini was equally di- 
vided between the central source and dust emission ac- 
cording to the estimated excess above the photosphere 
(CJ01). For the single-annulus models, the free param- 
eters of annulus width (2, 4, or 6 AU) and annulus in- 
ner edge distance from the star (2, 3, 4, 6, or 8 AU) 
are combined to yield fifteen models. We also generated 
models with two concentric and contiguous annuli with 
different but uniform brightnesses. In all of these mod- 
els, the inner annulus extends from 2 to 4 AU, whereas 
the outer annulus extends from 4 AU to either 6 AU or 
8 AU. The ratios of the total fluxes from the inner and 
outer annuli are either 1:1 or 3:1, respectively, resulting in 



four double-annuli models. The doublc-annulus is a sim- 
ple geometric construct that likely approximates a more 
complex radially-dependent brightness. The models were 
then convolved with a 2-D representation of the PSF that 
was generated from the azimuthally averaged PSF star 
observation at 18.3 /im. A ^-squared minimization com- 
parison of the resulting model profiles to the source data 
shows that the source is best approximated by a model of 
two annuli. In this model, a 2-4 AU annulus emits 75% 
of the thermal dust flux, and a 4-8 AU annulus emits the 
remaining 25%. Our models are not exhaustive but do 
illustrate the nature of satisfactory fits to the observed 
profiles. We also note that the dust emitting the unre- 
solved 10.4 /im excess must lie interior to the resolved 
18.3 /im flux and may be migrating inwards by P-R drag 
from the region of replenishment by collisions through a 
much lower-density region close to the star. 

3.2. Particle Size, Mass and Lifetimes 

The expected temperature for blackbody dust particles 
at 3.0 ±0.4 AU from the star is 320±H K. This is compa- 
rable to the color temperature of 3231^ K inferred from 
the two observed mid-infrared flux densities and associ- 
ated uncertainties. Thus, the particles may be blackbody 
emitters, in marked contrast to mid-infrared emitting 
grains in other debris disks; e.g., dust 50 AU from (3 Pic 
is at a temperature of >140 K compared to the black- 
body temperature of ~65 K (|Weinberger et al.1 120031 : 
iTelesco et aLll2005ft . Blackbody behavior requires parti- 
cle radii of a few microns or larger. Additional evidence 
that the particles orbiting £ Lep are large is provided by 
recent mid-infrared spectroscopy of the sourc e taken with 
the Spitzer Infrared Spectrograph (IRS) (|Chen et al.1 
2006). Silicates are commonly found in circumstellar 
disks, but the spectrum of £ Lep shows no evidence of 
a distinct silicate emission feature. Silicate particles as 
large as 1-2 /tm can result in a broad silicate feature, 
which is expected to disappear as particle radii increase 
above a few microns and the particles become optically 
thick in the mid-i nfrared (e.g., Przygodda et al. 2003). 
IChen et al.l (|2006h conclude that, if silicates are present, 
the majority of particle radii are larger than 10 /tm. 

We estimate the size of the smallest particles expected 
to remain on bound orbits {[3 — F r /F g < 0.5, Burns et 
al. 1979) in response to the forces of gravity and ra- 
diation pressure, assuming 14 Lq and 1.9 M Q for the 
stellar luminosity and mass, respectively, and unity for 
the radiation pressure efficiency. Adopting the standard 
silicate density value of 2.5 g cm~ 3 , we find that the min- 
imum radius is 3.4 /tm. Particles smaller than this are 
blown out of the system by radiation pressure within a 
few orbital timescales. The Keplerian orbital period at 
3 AU from £ Lep is ^3 years, so grains smaller than the 
blowout radius of 3.4 /im should be removed on a decade 
timescale. With this minimum particle size threshold, we 
estimate the number-density-weighted average dust par- 
ticle radius of 4.7 /tm by assuming that the incremental 
number of grains per radius interval da is proportional to 
a -3 5 , which describes the equilibrium size distribution 
of the dust part icle population in a collisional cascade 
(|Dohnanvilll969t ). 

Assuming that all particles have the average radius 
4.7 /tm and are in a thin shell 3.0 AU from the star, 
we estimate the dust mass from the amount of ultraviolet 
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and optical flux reprocessed as thermal emission (Jura et 
al. 1995, Eq. 5). With values for the fractional infrared 
lumin osity Lir/L* = 1.7 x 10~ 4 (|Faiardo-Acosta et al.1 
1998) and particle density 2.5 g cm -3 , we find a total 
dust mass of 6.7 x 10 21 g, or 0.37% of the total mass of 
the asteroids in our solar system. This value is ~40% of 
the mass estimate in CJ01 due to the smaller disk size 
(and therefore hotter dust) determined with our observa- 
tions. A change in the particle-size power-law exponent 
from 3.5 to 6.4 lowers the mass estimate by ~20%. 

Large dust particles (a > 3.4 /im) are su bject to loss 
by sp iraling into the star under P-R drag (jBurns et al.l 
1979). However, this process is only relevant if the 
time interval between destructive particl e collisions i s 
much longer than the P-R drag timescale (Wyatt 2005) . 
The characteristic collisional lifetime of the average ra- 
dius particle at a distance r from the star is t co \\(r) — 
iper(7")/(47TT G ff (r)), where t pcr is the Keplerian orbital 
period and T e g is the effective optical depth. We esti- 
mate the effective optical depth, defined as the ratio of 
the total cr oss-sectional area o f dust to the total area 
of the disk fWvatt et al.l Il999h . from the total area of 
average-sized (4.7 fj,m) particles necessary to produce the 
infrared luminosity and the area within the 2-8 AU an- 
nulus; for £ Lep, this value is 2 x 10~ 4 , approximately 
equal to the observed value for Lir/L*. Assuming a con- 
stant particle density, the ratio of the P-R drag timescale 
to the collisional timescale is a weak function of distance 
from the star, so we expect the relative significance of 
the dominant process to be similar throughout the disk. 
For example, the expected time for a particle to migrate 
by P-R drag from 6 AU to 3 AU (the center of each 
annulus) is ^16000 y, several times longer than the col- 
lisional timescale at 4 AU (the nominal division between 
the annuli) of only ~2000 years. Therefore, we consider 
it unlikely that the particles orbiting £ Lep migrate a 
substantial distance toward the star before being broken 
up in a collisional cascade. We conclude that initially 
radiation pressure and, subsequently, collisions are the 
two dominant processes defining the population and dis- 
tribution of particles in £ Lep's disk. 

4. DUST PRODUCTION SCENARIOS 

A plot of infrared excess versus age (Ric ke et al.ll2005l) 
of 266 A-type stars shows that ( Lep has a high infrared 
excess compared to other main-sequence st ars of simi- 
lar ag e. In this paper we have adopted the ISong et all 
(|2001f ) value of 231 Myr, estimated with photometry- 
determined stellar parameters to fit to stellar evolution- 
ary tracks; this estimate also incorporates the effects of 
rapid rotation, and a survey of A-type stars shows that 
£ Lep's position in the H-R diagram is close to that of 
other stars with comparable age estimates (|Jura et al.l 
119981) . However, if the age is 100 Myr or less (cf. Bar- 
rado y Navascues 1998; Song et al.), the excess level for 
the disk is within the range of values typical for those 
age bins. Therefore, the dust production scenario that 
we consider most probable depends, to some extent, on 
the age we assume for £ Lep: if the disk is older than 
100 Myr (cf. Lachaume et al. 1999; Song et al. 2001), 
the excess is indeed anomalously high and may indicate 
a recent cataclysmic collision, but, if the disk is much 
younger, we may be observing the dust production and 
decay associated with a relatively steady collision rate 



of smaller bodies. We examine both of these scenarios 
below. 

Recent cataclysmic collisions may acco unt for some 
of th e excess e mitting dust aro und /3 Pic (|Telesco et al.l 
2005) and Vega (|Su et al.l [2005). and may also contribute 
to the range of exc esses observed for stars of similar ages 
(|Rieke et al.ll2005l ). The destruction of one planetesimal 
of radius ~85 km could account for the entire observed 
mid-infrared-emitting dust mass in £ Lep. To assess the 
plausibility of a single cataclysmic collision of such a 
body resulting in a disk resembling our best-fit model, 
we consider the dynamics of the dust particles immedi- 
ately following their production in such an event. The 
initial orbital evolution of the newly formed dust par- 
ticles is dominated by radiation pressure, which acts to 
increase the semi-major axes of their orbits. If we assume 
that these particles are produced by the destruction of 
a body on a circular orbit at 2 AU (the inner edge of 
the inner annulus in the best-fit model to our observa- 
tions), we can estimate the radius for particles expected 
to migrate to a given distance, using standard relations 
(e.g., Kortenkamp & Dermott 1998; Wyatt et al. 1999). 
Since smaller particles will have larger apastron distances 
under the influence of radiation pressure, the size of a 
particle that can migrate to a given apastron distance 
corresponds to a lower limit on the size of particles in- 
terior to that location. The radius of particles expected 
to reach 4 AU (5.1 /im) will be both the upper limit for 
particles in the 4-8 AU annulus and the lower limit for 
the 2-4 AU annulus, where in the latter case all smaller 
grains have been blown farther out. The radius of parti- 
cles expected to reach 8 AU (3.9 /im) sets the lower limit 
for particle radii in the outer 4-8 AU annulus. 

To characterize the particle distribution that would 
yield the flux ratio (3:1) between the annuli of the model 
fit to our observations (§3.2), we sum over the areas of 
all the individual emitting particles, as described by a 
power-law particle distribution, assuming the radii are 
>5.1 /xm in the inner annulus and 3.9-5.1 /im in the 
outer annulus. (Since the radii of the particles that are 
expected to remain on bound orbits within 8 AU are 
larger than a few microns, we assume that the particles 
are blackbody emitters.) To estimate their temperature 
and flux, we then assume that all of the particles in the 
inner and outer annuli are at a characteristic distance of 
3 and 6 AU from the star, respectively. We find that, to 
match the model of our observations, the size distribu- 
tion of the dust population resulting from a collision at 
2 AU must be described by a p ower law w i th an exponent 
of 6.4, much steeper than the iDohnanvil (|1969fl value of 
3.5. Indeed, there is no reason to expect the immediate 
post-collision population to resemble the Dohnanyi dis- 
tribution that is thought to apply to a population that 
has evolved to collisional e quilibrium (|di Martino et al.l 
119901: iDurda fc Flvnrifl999f) . The post-collision distribu- 
tion presumably depends heavily on the details of the col- 
lision and the composition of the parent body. For exam- 
ple, with regard to the latter, "rubble piles" and differen- 
tiated asteroids may produ ce different relative amounts 
of small and large particles (|Grogan et al.|[200lh . A par- 
ticle size distribution consistent with our double- annulus 
model therefore suggests that the proportion of small 
particles produced in the hypothetical collision is larger 
than that implied by the Dohnanyi distribution. 
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In contrast to the catastrophic collision scenario, dust 
production by relatively "steady-state" grinding colli- 
sions may be more likely if the age of £ Lep is less than 
100 Myr and the excess level is not high compared to 
other stars of comparable age. Indeed, CJ01 propose 
that dust is replenished in £ Lep in an asteroid belt un- 
dergoing such a process. They estimate the total mass of 
the asteroids constituting the parent-bodies' mass by as- 
suming that the observed dust mass will, on average, be 
lost and replenished within the P-R lifetime of the par- 
ticles. If this has occurred continuously throughout the 
star's life, the required parent-body mass, Md x i a ge/ipR, 
is >4 x 10 26 g, which is ~5% of the primordial mass of 
asteroids (~ 1 M®, e.g., Lissauer 1987) in our solar sys- 
tem. However, since we contend that collisional cascades 
and subsequent radiation pressure effects are the domi- 
nant regulators of the dust population, we reconsider the 
estimate above for the more relevant collisional timescale 
t co u rather than that of P-R drag and determine a mass 
of ^6 x 10 26 g, ~10% of the solar system's initial aster- 
oidal mass. 

5. DISCUSSION 

There is a growing body of evidence for cataclysmic 
collisional events as a major, or even predominant, source 
of dust production in debris disks, as discussed in §4. 
A plot of simulated dust pr oduction rate in our solar 
system ([Dermott et al.ll2002l ) shows a slow decrease of 
dust area over time, punctuated by spikes from dust- 
producing collisions. This illustrates that the amount of 
dust alone cannot indicate how long it has been since 
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Fig. 1. — Profiles of azimuthally averaged normalized intensity for £ Lep (diamonds) and reference PSF star (dots); vertical lines indicate 
the FWHM values of profile fits to the PSF star and £ Lep. See the electronic edition of the Journal for a color version of this figure. 



